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To study the localization of random heteropolymers at an interface separat- 
ing two selective solvents within the model of Garel, Huse, Leibler and Orland, 
Europhys. Lett. 8 9 (1989), we propose a disorder-dependent real space renormal- 
. . . ization approach. This approach allows to recover that a chain with a symmetric 

0^ \ distribution in hydrophobic/hydrophilic components is localized at any tempera- 

Oi ■ ture in the thermodynamic limit, whereas a dissymmetric distribution in hydropho- 

bic/hydrophilic components leads to a delocalization phase transition. It yields in 
addition explicit expressions for thermodynamic quantities as well as a very detailed 
^ ■ description of the statistical properties of the heteropolymer conformations in the 

3 ' high temperature limit. In particular, scaling distributions are given for the lengths 

of the blobs in each solvent, for the polymer density, and for some correlation func- 
f*"- ' tions. In the case of a small dissymmetry in hydrophobic/hydrophilic components, 

the renormalization approach yields explicit expressions for the delocalization tran- 
sition temperature and for the critical behaviors of various quantities : in particular, 
the free energy presents an essential singularity at the transition (the transition is 
thus of infinite order), the typical length of blobs in the preferred solvent diverges 
with an essential singularity, whereas the typical length of blobs in the other sol- 
^.^ ' vent diverges algebraically. Finite-size properties are also characterized in details in 

fvj , both cases. In particular, we give the probability distribution of the delocalization 

f~^ ' temperature for the ensemble of random chains of finite (large) length L, and the 

^\ ' distribution of the numbers of blobs for the chains that are still localized at a given 

0^ ■ temperature. Finally, we discuss the non-equilibrium dynamics at temperature T 

starting from a zero-temperature initial condition. 
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Among the various subjects concerning the physics of polymers, the behavior of heteropolymers 
containing hydrophobic and hydrophilic components in solvents are of particular interest since they 
^ ' have obvious importance in biology ||l[. It is well known for instance that in a polar solvent, these 

heteropolymers prefer conformations where the hydrophilic components are in contact with the 
polar solvent, whereas hydrophobic components avoid contacts with the solvent. The behavior of 
heteropolymers in the presence of an interface separating two selective solvents, one favorable to 
the hydrophobic components and the other to the hydrophilic components, is less obvious, and 
has been much studied recently. In the pioneering work of Garel, Huse, Leibler and Orland M, a 
model was proposed and studied via Imry-Ma arguments, an analysis of the replica Hamiltonian 
and numerics : it was found that a chain with a symmetric distribution in hydrophobic/hydrophilic 
components is always localized around the interface at any temperature (in the thermodynamic 
limit), whereas a chain with a dissymmetric distribution in hydrophobic/hydrophilic components 
presents a phase transition separating a localized phase at low temperatures from a delocalized 
phase into the most favorable solvent at high temperatures. Experimentally, the presence of 
copolymers was found to stabilize the interface between the two immiscible solvents ||3|, since 
the localization of the heteropolymers at the interface reduces the surface tension. By now, the 
predictions of Reference |y] have been confirmed in the physics community by various approaches 
including molecular dynamics simulations Q, Monte Carlo studies 1^, variational methods for 
the replica Hamiltonian [^ 0], and exact bounds for the free-energy [||. Mathematicians have 
also been interested in this model, but the exact results obtained up to now are still far from a 
complete explicit solution. The localization at all temperatures for the symmetric case was proven 
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in References [|| jl^l by various criteria of localization on the paths measure. In the dissymmetric 
case, the existence of a transition line in temperature vs dissymmetry plane was proven in Reference 
PH l where some bounds for the free-energy were obtained in the asymptotic regimes of low and 
high temperatures. Finally, the relations between various localization criteria concerning either 
free-energy or path properties have been studied in Reference p^ . 

In this paper, we propose a new approach to study the model of Garel et al. M, based on a 
disorder-dependent real space renormalization procedure. Since this type of approach is not usual 
in the context of the thermodynamics of classical disordered systems, it seems useful to recall 
how these methods have been successfully applied in other domains. Disordered-dependent real 
space renormalization approaches have first been introduced in the field of disordered quantum 
spin chains ||I3] JMJ ||l^ p6| p7||. More recently they have been used to study random walkers 
in ID disordered environments |18| p^ , reaction-diffusion in ID disordered environments [Q and 
non-equilibrium dynamics of disordered classical spin chains [[l8[ [pi| . A similar method has also 
been introduced independently to study the coarsening of the pure one-dimensional $* model at 



zero temperature |22 . In all these fields, these methods, which are not "exact from first principles" , 
have however been remarkably successful in reproducing the exact results that were already known 
(see References Q Q ) , and in producing a lot of novel exact results for exponents and scaling 
functions for a large variety of physical quantities. These results apply in the large renormalization 
scale regime, corresponding respectively to low-tempcraturc in the field of disordered quantum spin 
chains, and to large time behaviors in the other fields concerning dynamics. These new predictions 
have moreover been confirmed numerically whenever they have been tested |15) ||2^ ||2J] . 

For the random heteropolymer problem, the disorder-dependent real space renormalization pro- 
cedure that we propose in this paper is defined to select the "important" configurations for the 
thermodynamics of a chain with a given realization of the disordered sequence in hydropho- 
bic/hydrophilic components at a given temperature. We expect this renormalization approach 
to become accurate in the large renormalization scale regime, corresponding here to the region of 
high temperatures. Since the localized phase extends up to a critical temperature Tc that goes to 
infinity as the dissymmetry parameter of the distribution in hydrophobic/hydrophilic components 
goes to zero, our approach allows to study the localized phases of symmetric chains and slightly 
dissymmetric chains. 

The paper is organized as follows. In the remainder of the Introduction, we recall the model 
of Garel et al. B in Section p^ A| , and give a summary of the main results of the present paper 
in Section IB. In Section O, we recall the physics of the heteropolymer problem at T = and 



at high temperatures, in order to motivate the introduction of a disorder-dependent real-space 



renormalization procedure in Section III. The results of the renormalization approach are given 
for symmetric and dissymmetric chains in Sections IV and [V| respectively. In Section [V|, we 



characteriz e th e finite-size properties of the problem by considering chains of finite (large) length. 



In Section VII , we discuss how the renormalization procedure also describes the non-equilibrium 
dynamics of a random chai n at high temperature starting at i = from a zero-temperature initial 
condition. Finally Section VII] contains the conclusions. In the Appendix [A[ we compare some 



results of the renormalization approach with the corresponding results obtained by Garel and 
Orland via an extension of the usual Imry-Ma argument 



A. Model and notations 

We consider the model introduced by Garel et al. g with some slightly different notations: a 
polymer chain consists of L monomers indexed by i — 1,2, . . . , L, where each monomer i carries 
a quenched random charge qi. The {qi} are independent identical random variables drawn with 
some probability distribution C{q). An interface at z = separates a solvent in the domain z > 
favorable to positive charges q > 0, from a solvent in the domain z < favorable to negative charges 
g < 0. More precisely, the energy of a configuration C — {fi = {xi,yi, Zi)} of the heteropolymer 
characterized by a realization {gi, (72, • ■ ■ ^l} reads 

L 

E{C) ^ -Y,q,sgn{zi) , (1.1) 

i=l 



and the partition function reads 

L-l 



ZL{P;{qi}) = Tra,ce[p.yY{5{\ri+i-ri\- a)expip^qiSgn{zi)\ . (1.2) 



1=1 



In a continuum version of this model, the nionomer index i becomes a continuous variable s, and 
the {x, y) coordinates play no role so that the partition function becomes a path integral over 
one-dimensional Brownian trajectories {z(s)} [^1 



Zl{P-A<i{s)}) 



Jvzis)e^pL^£ds(^^^ + (3 £ dsqis)sgn{zis))\ . (1.3) 



In this paper, we propose a disorder-dependent real-space renormalization procedure, that should 
give valid results in the universal regime of large scales, where the details of the microscopic charge 
distribution C{q) are not important, as long as the distribution C{q) satisfy the hypothesis of the 
Central Limit Theorem. The important parameters of C{q) will thus be the mean value and the 
variance 

dq q C{q) = qo (1.4) 

-OC 

/ + 00 
dq{q-qo)^ C{q)^2a . (1.5) 

-oo 

To study the biased case qo j^ 0, we will always choose the convention qo > and introduce the 
parameter S defined as the non-vanishing root of the equation [|l9| 

/ + OC- 
dqe-2*?C(g) = 1 . (1.6) 

-OO 

For instance, in the case of a Gaussian distribution considered in References [@1 @ [|| 

C{q) = -l^e^''^ , (1.7) 

the parameter S is simply the ratio between the mean value go ^nd the variance (2ct) |14| 

..|i . (1.8) 

For the case of binary distribution considered in References Q ||ll|] |1J] where the charge q 
takes a positive value q-^. with probability c and a negative value (— ?-) with probability (1 — c) 
then 



go = cg+ - (1 - c)g_ , (1.9) 

cr="-^^{q+ + q-? , (1.10) 

and S is the solution of 



ce 



-2Sq 



+ + (1 - c)e+2''«" = 1 . (1.11) 



For instance in the case g_ = g-|_ = gi, we have 

go-(2c-l)gi (1.12) 

a = 2c{l-c)ql (1.13) 

6=±\n-^ . (1.14) 

2gi 1 — c 

An important property of the parameter 6 is that it becomes universal and coincides with the 
Gaussian expression in the limit of small dissymmetry (5^0 |l4| 



S^^ . (1.15) 

Concerning the spatial behavior of the polymer, the only important parameter will be the diffu- 
sion coefficient D characterizing the large distance behavior of the free chain 

<{z{s)~z{s')f >^D\s-s'\ . (1.16) 

We stress that in this paper, we do not take into account excluded volume effects because they do 
not have any crucial effect for the problem of the localization of the heteropolymer at the interface 

iii- 

B. Quantities of interest and summary of main results 

The quantities of interest to characterize the localized phase of heteropolymers around the 
interface are on one hand the free-energy f{T) per monomer that governs the thermodynamic, 
and on the other hand the statistical properties of the spatial conformations of the chain. In 
particular, the important characteristic lengths are the typical lengths luobO^) ^^ blobs in the 
(±) solvents as well as the typical distances z^i^^^iT) between the chain and the interface. At a 
more refined level, one may also consider the probability distributions P^{1) of the lengths of (±) 
blobs, the density p^{z) of polymer at a distance z from the interface in the (±) solvents, and 
even correlation functions as for instance the solvent-solvent correlation < sgn(2;(s))sgn(z(s')) >. 
Finally, we will also be interested in the finite size effects of the problem. 

Before introducing our approach and proceeding with the calculations of these quantities of in- 
terest, we now summarize the main results that will be derived from the real-space renormalization 
approach in the remainder of the paper. 



1. Summary of main results for the symmetric case (qo = 0) 

For the case of symmetric chains (qa = 0), the renormalization approach yields the following 
behaviors at high temperature 

hiobiT) {\nTf , (1.18) 
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z{T)^^{\nT) , (1.19) 



in agreement with the Imry-Ma argument results [0| (see also Section II B for the discussion of the 



presence of the logarithmic correction). In addition, we obtain that the rescaled length of blobs 
is distributed with the law 

OO / 1 \ 

P(A)= y («+- )7r(-l)"e-'^'^("+5)' ~ Tre-'r^ (1.21) 

n— — oo ^ ' 

oo 

^ T7l — — oo ^ 

The polymer density takes the scaling form p{z)dz = R{Z)dZ where the rescaled distance to the 
interface 



is distributed with the following law R{Z) that decays exponentially at large distance 

R{Z) ^4 dXP{X)VX / dwe~"' ~ -V2Ze~''^ . (1.24) 

Jo J-S= ^^°° "■ 



n/T 



Finally, to characterize the finite size effects, we will compute the distribution of the delocalization 
temperature Tdcioc over the ensemble of random cyclic finite chains of (large) length L. The result 
is that the rescaled variable 
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is distributed with the law 

2 +°° 

^ n=l 
o 4-00 



2 — ■> " ^ 

71 e ■*9 ~ — r-e ^is 



-^ y [2(2m+l)V-l]e-('"'+i)^ ^ -^ 

* m ^ — rv-i * 



o ^ 

5'e-^ 



(1.25) 

(1.26) 

(1.27) 



771 — — OO 



^. Summary of m,ain results for the dissym,metric case (go > Q) 



For the case of dissymmetric chains (go > 0), in the limit a '^ qq, the renormalization approach 
yields that the delocalization transition takes place at the critical temperature 



T, 



4cr 
3qo 



(1.28) 



in agreement with the scaling obtained previously by other methods [0 B H] . In addition, we 
obtain that the transition is of infinite order, with the following essential singularity for the free 
energy the free-energy /(T) 



/(T)-/(T,) ~ -2gofln— )cxp 



In^ 

90 



(1.29) 



For the statistical properties of the heteropolymer chain, we find that the typical length ?^o(,(r) of 
blobs in the preferred solvent diverges with an essential singularity at the transition, whereas the 
typical length l^i^i,{T) in the other solvent diverges algebraically 



^tiob^T) ^ -2exp 



In^ 
go 



luobiT) 



In- 



T^T- 95 ( 1 - ^ 



The rescaled length of blobs in the preferred solvent, defined as 



T^exp 

'Jo 



go 



(1.30) 
(1.31) 

(1.32) 



is distributed with the exponential law e ^+ near the transition. The polymer density p'^{z) takes 
the scaling form p^{z)dz = R^{Z)dZ where the rescaled distance to the interface 

(1.33) 



2 90 



exp 



In; 



2(1-^) 



is distributed with the following scaling distribution that decays exponentially at large distance 



R+{Z)^A due"'' / dvv^e-'' ~ VvrZe^^^ . (1.34) 



n ./A ^ — *oo 



Finally, to characterize the finite size effects, we will compute the distribution of the delocalization 
temperature Tdcioc over the ensemble of random cyclic finite chains of (large) length L. The result 
is that the random variable 



a / 4(7 \ ^c-^doio 



r 



IqL V ii 



(1.35) 



is distributed with the law 



D+{r) = —e-r . (1.36) 

In particular, the typical value for the delocalization temperature presents a correction of order 
(1/lnL) with respect the critical temperature Tc 



II. PHYSICAL PICTURE AT T = AND AT HIGH TEMPERATURES 

In this section, we recall the physics of the heteropolymer problem at T = and the Imry-Ma 
picture of Reference valid at high temperature. They will be useful to motivate the introduction 
of the real-space renormalization procedure in the next section. 

A. Description of Ground States 

At T ~ 0, for a given realization {qi, q2, . . .ql}, each monomer i will be in its preferred solvent 
according to 

sgn(zj) = sgn(qj) . (2.1) 

The ground states have thus for energy 

L 

Eo = -Y,W^ . (2.2) 

i=l 

The corresponding configurations of the chain are all the random walks that cross the interface each 
time there is a change of signs in the realization {gi, q2, ■ ■ ■ Ql}- More precisely, it is convenient to 
decompose the reahzation {gi, (72, ■ • • 9l} into groups (a) of consecutive charges of the same sign, 
each group (a) containing a number l^ — 1, 2, ... of monomers, and carrying an absolute charge 
Qa with 

1<. = Y.^ (2-3) 

Qa = 1^1*1 • (2-4) 



Then the ground states of the chain can be decomposed into blobs a containing I a ^ 0(1) 
monomers. The chain is thus localized around the interface with a typical distance of order 0(1). 



B. Imry-Ma arguments at high temperature [2] 



In this section, we recall the Imry-Ma arguments |26| of Garel et al. g , for the symmetric and 



dissymmetric cases respectively. They are expected to be valid at high temperature. 



1. Imry-Ma argument for the symmetric case 

The Imry-Ma argument exposed in for the symmetric case (go — 0) can be summarized as 
follows. Assuming that the chain is localized around the interface with typical blobs of length I in 
each solvent, the typical energy gain per blob scales as -s/ol, whereas the reduction of entropy per 
blob scales as InZ (in H, only powers of I were considered for the symmetric case, and thus In^ was 
replaced by Z° ^ 1, but since this \nl dependence plays a crucial role in the dissymmetric case, we 
prefer to keep this In I dependence everywhere in this paper) . Neglecting prefactors of order 1 , the 
optimization of the free-energy per monomer 



with respect to I leads to 



I r2 



(InO^ a 



(2.6) 



At high temperature, the scaling of the typical blob length 1{T) and of the free-energy /(T) are 
therefore given by : 

/m^-j^ (2.8) 

and thus the heteropolymer remains localized at any finite temperature. 



2. Imry-Ma argument for the dissymmetric case 

The Imry-Ma argument exposed in for the dissymmetric case is actually much more sub- 
tle than in the symmetric case, because to describe correctly the blobs in the (— ) solvent, it is 
necessary to consider the "rare events" where the sum of random variables qi of positive mean 
gi" = go > happens to be negative enough in order to make more favorable for the heteropoly- 
mer to make an excursion in the (— ) solvent rather than to stay in the otherwise preferred (-I-) 
solvent. More precisely, the Imry-Ma argument of Reference g is as follows in our notations : at 
high temperature, the heteropolymer is expected to be mostly in the preferred (-I-) solvent, except 
when a blob of length l^ in the (~) solvent becomes energetically favorable, with a blob energy 

Q- = —J2i=j 94 > that is "sufficient". The probability to have X]i=j 1i — ~Q~ is given by 
the Gaussian 

1 (O'+qQi-)" 

Prob(Q-) ^ e T^ . (2.9) 

Thus the typical spacing l^ between two such events behaves as the inverse of this small probability 

l+^e J^F . (2.10) 

As a consequence, in a configuration with blobs of order (/+,^^) with ^+ ^ ^^, the energy gain 
Q~ due to one excursion in the (— ) solvent behaves as 



^/ Aal- \nl+ - q^r , (2.11) 



whereas the loss of entropy due to this excursion is of order lnZ+ + ln/~ ~ \nl~^. The free-energy 
difference per monomer between this locahzed state and the delocahzed state in the preferred 
solvent can then be estimated as (neglecting prefactors of order 1) : 

1 . „_ ^, ,^x 1 



f{T,l+,l-) - fdeiociT) ~ ^ (-Q- +rin;+) ~ — i^qol- - V4al- \nl+ + Tlnl+j . (2.12) 

The optimization with respect to l~ leads to 



(T 



r ^ —\nl+ (2.13) 

% 



and 



Q- ln/+ (2.14) 

Thus both the energy gain Q~ and the entropy cost have the same ln/+ dependence! As a 
consequence, the free-energy difference factorizes into 

/(r,/+)-/rfe/oc(r)^(r-T,)i^ , (2.15) 



where the parameter 



Tc^- (2.16) 

Qo 



thus represents the critical temperature between the localized phase T < T^ and the delocalized 
phase T > T(.. Contrary to the symmetric case, the behavior of the free-energy with respect to 
the temperature remains unknown, since the typical blob length l'^ depends upon the temperature 



T in a way that is not determined by the present argument. Still, the relations (2.13) and (2.15) 
between the typical blobs lengths l'^ and l~ and the free-energy f{T) represent non-trivial tests 
for any theory constructed to describe the heteropolymer problem. 



III. DEFINITION OF AN EFFECTIVE THERMODYNAMICS 

To go beyond the Imry-Ma arguments of previous section, one needs to consider probability 
distributions of blob lengths and blob energies. In Appendix ^, we reproduce an extension of 
the usual Imry-Ma argument which was introduced by Garel and Orland [ p5[ , and which shows 
how to get the asymptotic behavior of the probability distribution P{1) of the blob length I in the 
limit of small I. Here, to study the probability distributions of blob lengths and blob energies, we 
will construct via a real-space renormalization procedure the "optimal" Imry-Ma domain structure 
associated with a given heteropolymer at a given temperature. We will then argue that the 
configurations corresponding to this optimal Imry-Ma domain structure dominate in the partition 
function asymptotically at high temperature. 



Definition of a Real-Space Renormalization procedure to construct the optimal 

Imry-Ma domain structure 



At T = 0, we have seen in Section [I A that the chain is decomposed into blobs a containing 
la consecutive monomers of charges of the same sign, and carrying an absolute charge Qa ■ As T 
grows from T = 0, we consider the configurations of the chain obtained from the ground states 
structure by iteratively flipping the blobs of smallest absolute charge Qmin = T- When we flip the 
blob {Q2 = r, I2) surrounded by the two neighboring blobs (Qi, li) and (Qa, Z3), we obtain a new 
blob of absolute charge Q and length / given by (see Figure 



(^iM) 



Q^Qi + Q3-Q2 

l = ll+l2 + h ■ (3.1) 



(Q,,l,) 



(Qyl,) 



z=Q 




(Q,,l,) 



iy 






z=0 




FIG. 1. Illustration of the real-space renormalization procedure : the flipping of the blob (Q2 ~ XjI?) 
surrounded by the two neighboring blobs {Qi,li) and (Qsja) gives a new blob {Q,l) with the rules ([3.l|). 



The renormalization procedure corresponding to the rules (3.1) has already been extensively 
studied in the context of the Random Transverse Field Ising Chain |lj|, where the role of the 
absolute charges Qi was played by the logarithm of the random couplings and random fields (In Ji 
and In/ii), and in the context of Random walks in disordered ID environments []l9| where the 
role of the absolute charges Qi was played by the energy barriers Fi. As explained in details 
in Reference |19|] , the renormalization procedure constructs iteratively the large scale extrema 
statistics of random walks, where only barriers bigger than a given scale F are kept. Here the 
corresponding random walk is simply given by the sum of the quenched charges ^q Qj as a function 
of the monomer index i, and the renormalization procedure gives the "best" blob structure, given 
the constraint that only blobs of charges bigger than a given scale F are kept. 

To establish a relation between the renormalization scale F and the temperature T, we now 
have to determine the conditions under which the flip of the blob {Q2, h) is favorable. The cost in 
energy of the flip of the blob {Q2J2) is simply 



whereas the corresponding gain in entropy reads 

AS^^'P = ln(M(/i +I2 + h)) - \n[Mili)M{l2)M{l3)] 



(3.2) 



(3.3) 



where A4{1) represents the number of ID random walks of I steps going from z = to z = in 
the presence of an absorbing boundary at z = 0^. As a consequence, the free-energy difference 
corresponding to this flip reads 



/^pfl^P = /s^Efl^P _ TASf^'P = 2Q2 - Tin 



/ Mjh+h + h) 
\Mih)M{l2)Mih) 



(3.4) 



The optimal Imry-Ma domain structure corresponding to a given realization of the disorder 
at a given temperature is thus constructed as follows : we iteratively flip the blobs of smallest 
absolute charge Q, starting from the ground-states blob structure, and continue as long as these 
flippings produce a decrease of the free-energy [AF^'-^p < 0). We stop the procedure when the 
next flipping would correspond to a raise of free-energy [AF^^^p > 0). We will define Feq(T) as the 
renormalization scale F where we have to stop the renormalization. We now present a dynamical 
interpretation of the renormalization procedure. 



B. Dynamical interpretation of the renormalization procedure 

The renormalization procedure defined above to construct the optimal blob structure at tem- 
perature T can be given a perhaps more direct physical meaning if one considers the dynamics at 
high temperature T for i > starting from a zero-temperature initial condition at t = 0, i.e. after 
a quench to T = for t < 0. Indeed, the time necessary to flip a blob of absolute charge Q follows 
an Arrhenius law t ^ toe^^^ (where to is a microscopic time scale). The dynamics thus corresponds 
to the iterative flipping of the blobs of smallest absolute charge remaining in the chain described 
by the renormalization procedure, where the renormalization scale F now corresponds to time via 

r = Tln— (3.5) 

to 

as in References ]l8| |l9[ . The dynamics takes place up to time teg where equilibrium at temperature 
T is reached 

Tin^ = re,(r) . (3.6) 

to 



C. Discussion of the vaHdity of the effective thermodynamics 

Let us deflne the effective partition function Z£*(/3;{qi}) as the sum over all the configura- 
tions that correspond to the disorder and temperature dependent optimal Imry-Ma domain struc- 
ture constructed via the renormalization procedure defined above. The real partition function 
Zl{(3; {qi}) can be decomposed into 

Zl{P; {qi}) — Z'j^{f3; {qi}) + sum corresponding to other Imry-Ma domains structures . (3.7) 

The effective thermodynamics based on Z'^{f3;{qi}) is of course very approximate at low tem- 
peratures. However, at high temperatures, the renormalization scale F = Qmin will be large and 
the probability distributions of absolute charges of the blobs will become infinitely broad WM. 
More precisely in the symmetric case, the difference between the typical value Qtyp of the absolute 
charges and the minimum value Qmin — T will also be large and of order F : Qtyp — Qmin ~ T- 
As a consequence, we expect that all the configurations that do not correspond approximativcly 
to the optimal Imry-Ma domain structure will be highly suppressed in the partition function, in 
comparison with the "optimal configurations". The validity of the effective thermodynamics is 
thus based on the idea that, for a given heteropolymer at a given high temperature, the large scale 
Imry-Ma domain structure is essentially unique. 

As explained in the Introduction, in the other fields where disordered-dependent real space 
renormalization approaches have been used, they have been able to reproduce the non-trivial 
exact results obtained previously via other methods (see References ||lj] [|l^ ). These agreements 
with the Mc Coy-Wu exact results in the field of disordered quantum spin chains [|^ , and with the 
Kesten-Golosov distribution in the field of ID random walks in random media p9| | have provided the 
best evidence for the asymptotic exactness of the renormalization approach in these fields, and give 
confidence in all its other novel predictions. Here for the heteropolymer problem, we do not have at 
our disposal so precise exact analytical results derived via other methods to "test" unambiguously 
the asymptotic exactness of the predictions of the renormalization approach. As a consequence, we 
hope that in the future, the explicit results of the renormalization approach presented in this paper 
will be tested precisely by numerical studies. It would be also very interesting to test more directly 
the idea of the dominance of the optimal Imry-Ma domain structure in the partition function, by 
comparing, sample by sample, the typical configurations of the heteropolymer at equilibrium at 
temperature T with the optimal blob structure obtained by the numerical implementation of the 
renormalization procedure up to scale Feg(r). 

The remainder of the paper is devoted to the detailed study of the properties of the effective 
thermodynamics. To proceed now with the concrete calculations of the renormalization approach, 
it is again more convenient to discuss separately the symmetric case <7o = smd the dissymmetric 
case qo > 0. 
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IV. STUDY OF THE SYMMETRIC CASE Qo = 

A. Properties of the renormalization procedure 

It was shown in Reference jlj] that with the above renormahzation procedure (^J), the proba- 
bility distribution Pr{Q,l) that a blob at renormalization scale F has absolute charge Q and length 
I flows towards a fixed point distribution P*{ri, A) in the rescaled variables 

V^^^^jr- and X ^ -^ , (4.1) 

with a defined in ( |l.4| ), and that this fixed point distribution P*{rj, A) reads in Laplace transform 
with respect to A 

£a^,(P*(77,A)) = r dXe-'^P*{rj,X) - -^^ g-.^icothyi _ (43) 

Jo smhVs 

In particular, the distribution of the rescaled absolute charge 77 is a simple exponential 

/■OO 

P*{V)= / rfAF*(r;,A)=e-" , (4.3) 

Jo 

and the distribution of the rescaled length A of a blob reads 



f°° / 1 \ f^''°° ds e"^ 

n— — 00 ^ ' ^ r?i— — 00 

We also introduce the notation Pflip(A) for the distribution of the bonds which are about to be 
decimated ( i.e. having ?7 = ) 



(4.4) 



-Fflip(A) — £ ;i^ 1 . , _, 

\smh(vs) 



•n — — rvo * ryyi — — 1 



X(™+2)'-^ 



e"i("+^)' . (4.5) 



B. Relation between temperature and renormalization scale 

We now examine the condition APf'-^P < 0, where the free-energy difference AF^^'^p due to a 



blob fiip was given in Equation (3.4). For large I, the number M{1) of ID random walks of I steps 



going from z = to z = in the presence of an absorbing boundary at z = behaves as 

M{l)^K-j^ , (4.6) 

where k is a constant, and where In/i characterizes the entropy of the free random walk (for 
instance /j, = 2 for the ID lattice). In the effective thermodynamics at large renormalization scale 
F, the free-energy difference of a blob flip then reads 

AF/'^^ = 2g,-rinf— M^^l!^) . (4.7) 

yn^ih + h + hy/^j 

Using the rescaled variables A^ = y^, we obtain the average 



AF/'^P = 2r-3Tln— +T(-X + lnK2) , (4.8) 

a 2 
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where _ftr is a pure numerical constant that can be obtained from the distributions (4.4) and (4.5) 



K= 



dXiP*iX,) / d\2P^.(\2) / dA3P*(A3)hl(Ai + A2 + A3) 



/•OO /'CX3 

-2/ dXP*{X)\nX- dAPfl%(A)lnA 

Jo Jo 



(4.9) 



The scale Teq{T) where the renormalization procedure has to be stopped is the solution of 
AF^^^P(r) = and drAF^^^P{r) > 0. For large F, and thus for large temperature, we find that 
Feq(r) is implicitly defined as the inverse of the function r(Feg) given by 



nve,) 



eq 



31nF 



(4.10) 



eq 



As a consequence, for any arbitrary large temperature T, we find that the renormalization 
procedure has to be stopped at the finite large renormahzation scale Feg(r). The polymer is 
thus always localized around the interface, with blobs of typical length l^^°^[T) behaving at high 
temperature as 



hiob{T) 



F 



q T^{\nTf 



a a 

In this regime, the typical distance z{T) to the interface behaves as 



(4.11) 



D, 



z{T)^^DlMob{T)^X -T. 



eq 



'£.T(\llT) 

a 



(4.12) 



The se scaling behaviors are in agreement with the Imry-Ma argument of Reference Q (see Section 
II B , where the logarithmic correction is found to be present in the Imry-Ma argument) and with 
the exact free energy bounds of Reference [g| . The Replica variational approaches of References 
Q give the same scaling behaviors without the logarithmic correction. 

We now study more precisely various physical quantities for which the present renormalization 
approach yields explicit predictions. 



C. Thermodynamic quantities 



The energy of a blob is simply given by minus its absolute charge (— Q), and thus the energy per 
monomer e(T) of the chain in the thermodynamic limit can be obtained through a decomposition 
into blobs indexed by a as 



s(T) = 



EaQ«_ /Q,/Q^r.,(T)(Q,0 



l^a '■o 



InJPr.,iT){Q,l) 



(4.13) 



with Feq(T) defined in (4.10). Using the fixed point solution of Equations (4.1, 4.2), we obtain at 
large temperature 



,{T) c^ - 



4cr 



4ct 



re,(T) 



STlnT 



(4.14) 



The entropy associated with the z direction for a blob of length / is given by ln(A^(/)), where 
■A4{1) behaves as in Equation (4.6) in the large I regime. As a consequence, the entropy per 
monomer s{T) is given in the thermodynamic limit by 



s(T) = 



EJ^alnM + lnK-|lny _ /q JHn/x + In^ - | lnZ]Pr^^(j,)(Q, Z) 



l^a '■a 



ki^Pr.^miQJ) 



(4.15) 



Using again the fixed point solution (4.1 , 4.2), we have at large temperature 



12 



s{T)c:i In/i 



GalniVeqiT)) ^ a 



r?,(r) Ti^{T) 



/•OO 

31no- + 21nK-3 / dXP*{X)lnX 



Finally, we obtain the free-energy per monomer /(T) of the chain at large temperature 

The dominant behavior depending on the disorder parameter a is in agreement with the Imry-Ma 
argument of Reference [EJ as explained in Section [IB. 



D. Distribution of the blob lengths 

In the renormahzation approach, the lengths of the different blobs are independent identical 
random variables, and the rescaled blob length defined as 

A^^ (418, 

is asymptotically distributed at high temperature with the distribution P*(A) given in Equation 



(4.4). In particular, its asymptotic behaviors for large and small A are respectively given by 



P(A) ~ Tre""^^ (4.19) 



''^'K-ov^^'^ ■ ^'-''^ 



As noted by Garel and Orland [g5|, the small A behavior found in Equation (4.20) for the 
distribution P{X) can be obtained via an extension of the usual Imry-Ma argument, as explained 
in the Appendix (H). 

E. Density p{z) of polymer at a distance z from the interface 

We now study the density p{z) of polymer at a distance z from the interface as follows. We first 
need to introduce the probability St {I, x) that a given point of the chain belongs at renormahzation 
scale r to a blob of length / and is at distances (a;, / — a;) from the ends of the blob. It reads 

Sril,x)= J\^\\ e{x)e{l ~ x) , (4.21) 

Jo dl I Pt{1) 

and is normalized to /p dl /g dxSr{l, x) = I. We then need the probability Bx.i{z) to be at height 
z at "time" x for a Brownian motion of diffusion coefficient D starting at z = at x = and 
finishing at z = at a; = Hn the presence of a reflexive boundary at z = 

2 / ; xV2 



It is normalized to /„ dzBx.i{z) — 1. The probability p{z) that a given point of the polymer is 
at a distance z from the interface (normalized to L dzp{z) — 1) can be now expressed within the 
renormahzation picture as 
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/•OO pi 

p{z)= / dl dxSr^^{T){l,x)Bx^i{z) 

JQ Jo 

Introducing again the rescaled variable A = ^ , we thus obtain the scahng form 

p{z)dz ^ R{Z)dZ , (4.24) 

where the rescaled distance Z to the interface, defined as 

r2a z 



D re,(T) 

is distributed with the following scaling function 



(4.25) 



poo /'OO 

R{Z)^aI d\P*{\)\fX I due-""^ . (4.26) 

where P*(A) is the fixed point distribution of Equation ( |4.4| ). Deriving with respect to Z yields 

R'{Z) = -4 / rfAP*(A)e"^ = -16Z^(-l)"ifi [{2n + 1)ttZ] (4.27) 

in terms of the Bessel function Ki (y) . 

In particular, the asymptotic behaviors at Z — > and Z —^ oo are given respectively by 

i?(Z) =i?(0)~4Z + O(Z2) , (4.28) 

RiO)^2V^l dAP*(A)VA.-g^i=^ , (4.29) 

and 

R{Z) ~ -%/2Ze-'^^ . (4.30) 

Z — ^oo 7T 

We may also compute the moments 

r(i + |) '■°° 



/ dZZ''R(Z)=2^—^ dXP*(X)X^+i 
Jo k + 1 Jq 



,r(i + |)r(2 + |) - (-1) - 

n— ^ ^ ■' 

111 particular, the mean value of the square distance to the interface is given at high temperature 



(fc + l)7r3+fe ^ („+ i)3+fe • ^ -^ ^ 



by 

-^z^ ^ rl,{T)- j^ dZZ^RiZ) = ^1,{T) ^ —T^lnTf . (4.32) 

F. Density pa.biz,z') for two thermal copies of the same chain 

To compare our approach with the Replica Gaussian Variational description of Reference H], 
we now compute the probability pa.b{z,z') that a given monomer of a given polymer chain is at 
distance z from the interface in a configuration (a) and at a distance z' from the interface in a 
configuration (b) (note that in the effective thermodynamics, we have sgn(z) = sgn(z')). Using 
the same notations as in the previous section, the joint distribution of z and z' can be obtained as 

14 



f OO pi 

Pa,b{z,z')= I dl dxSr^^(T){l,x)Bj;j{z)Bx^i{z') 





Using again the rescaled variable A = ^ , we obtain the scahng form 

PaAz,z')dzdz' ^Ra,biZ,Z')dZdZ' with ^ = J ^ —^ . (4.34) 

The scaling function Ra,b{Z^ Z') reads in terms of the fixed point distribution P*(A) of Equation 



(4.4) 



RaAZ;Z') = -j dXP*{X)e i^Ko i ^^ ) . (4.35) 

It is interesting to note that it is actuaUy a function of the single variable (Z^ + Z'^). 

The probability that a given monomer is at the same z-coordinate in the two configurations (a) 
and (b) thus behaves at high temperature as 

/ dzpa.b{z,z) = J — ——— dZRa.b{Z,Z)^ — — — , (4.36) 

as compared to the decay as (l/T) found in Reference jj] for the Replica symmetric solution. 



G. Solvent-solvent correlation function < sgn(z(s))sgn(2(s')) > 



Within the renormalization picture, the correlation function < sgn(z(s))sgn(2:(s')) > of the sol- 
vents seen by two monomers at distance (s — s') along the polymer chain turns out to correspond 
exactly to the spin correlation function {SQ(t)Sx{t)) of the random field Ising model in the Glauber 
dynamics starting from a random initial condition computed in References |18| pi| . We thus only 
quote the result 



< sgn(z(s))sgn(z(s')) >= C ^ ^,_,,^ 



- uaiiii I 
P P 

.\s-s'\ 



— tanh 



^ (2n+ 1)4^4 e r , ^4.d/j 

n— — C30 ^ ^ 

where F = TeqiT). In particular, the correlation length ^(T) associated with this correlation 
function is given by the decay of the dominant exponential at large separation \s — s'\ 

V. STUDY THE BIASED CASE Qo > 

A. Properties of the renormalization procedure 

In the biased case, it is necessary to introduce two probability distributions P^{Q,l) and 
Pf{Q,l) representing the probabilities for a blob at scale F in the (±) solvent to have absolute 
charge Q and length /. It was shown in Reference [H that with the renormalization procedure 



(3.1), these probability distributions flow at large F towards 
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{Q,l) - C'^li ([/±(p)e-W-r)"r*(rt) (5.1) 

u^{p) = A{p) coth [rA(p)] T S 



r ^^^ - sinh [rA(p)] 



Aip) = JS-- + P 



where a and 5 have been defi ned in Equations (L4) and (L£) in terms of the initial charge dis- 



tribution C{q). The solution s (5.1 ) are valid in the scaling regime of large F, small 5 and small p 
with 5T fixed and pT^ fixed |l4|. 



In particular, the distributions of the absolute charge Q >T are simple exponentials : 

25 



'+(g) = u+e-(Q-n< 


with 


Mp 


'r"(Q)-%e-('?-i^)"r 


with 


Wp 



(5.2) 



-2r5 ' (^•^) 



1-e 
whereas the distributions of the length of blobs have the form of an infinite series of exponentials 

El 

P±a)=c-' ( ^(^)^^^" 

r I ) '-p^i y^^p^ ^Qgj^ [rA(p)] T -5 sinh [rA(p)] 

OO 

-^5:j±(7 = ^r)e-^^J(-^^) , (5.4) 



r2 



n=0 



where the functions 5.^(7) and J^il) a-re defined in terms of the roots 0^(7) {n = 0, 1, ...) of the 
equation 

a^(7)cotan(a^(7)) = ±7 with nvr < a,'^(7) < (n + l)7r . (5.5) 

For 7 > 1, the root ckq (7) does not exist, but is replaced by the positive root fig (7) of the equation 
0^0^(7) coth((5(j"(7)) = 7. In terms of these roots, we have 

s^(7)-7' + K=^(7)f , (5.6) 

2(-l)"(a±(7))V^' + (""(7))'^^^ . . 

Jn (7) = .. L ... , (5.7) 

7 + (an(7)) T7 

except for the (+) n = term in the domain 7 > 1 for which 

4(7 > 1) - 7' - i&Ul)f (5.8) 

, 2(5„+(7))V^'-("^(^))'e"' , . 

•^o+(7 > 1) = ZT^Tr^2 2 ■ (5-9) 

7+(ao(7))^-7^ 

The mean lengths of the blobs in the domains (±) have simple expressions 



— J_ 

4(52 cr 



^r=77^(e''^-2<5F-l) , (5.10) 

— 1 



^^ = w^^^^^~^ + ^'''^^ ' (^-^^^ 



but the mean values of In Z* are unfortunately much more complicated : 



lnl±^ rdl{lnl)P^il)^f]:^^\n(^—)+\n(—)-CEuier , (5.12) 



where CeuIci is the Euler constant. 
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B. Relation between temperature and renormalization scale 



Using again Equation (|4.6|), the averaged free-energy difference corresponding to the flip of a 
blob out of the domain (±) reads 



/^pfHp{±) = 2r - Tin 



mtifr/' 



{lf+lf+lf)V2^2 



(5.13) 



which yields to a quite complicated expression for arbitrary j ^ 6T. 

In the following, we consider the case of large 7^1, where the distribution P^ is dominated 
by the 71 = term [111 



P+(0=ia+(7)e-"''('')' 



with 



a+(7)~4,5Ve-2T 



(5.14) 



leading to 



ln/^ 



In 



a+(7) 



CeuIci- = 27 + In 



Ad^c 



C, 



Eulor 



(5.15) 



In this regim e, th e probabilities to have blobs with absolute charge Q = F in the (±) domain read 
respectively (5.2) 



P+(g = r) =u+~2(5, 
p-{Q = r) = Up -2(5 



-2r6 



(5.16) 
(5.17) 



and thus the renormalization consists asymptotically in flipping only blobs out of the (— ) domain. 
We thus now concentrate on the free-energy difference corresponding to the flip of a (— ) blob, 
which can be estimated as follows 



AF/'»p(-)= 202-rin 



ii+i^itr^' 



(^l++l-+l+)V2^2 



3 

2r- -T 

2 

2r--T 

2 



21n(;+)-ln(?+ + /+) 



2Sr + In 



4(52 c 



-Q 



Eulcr 



- 1 



(5.18) 



so that finally, in the regime (5 — )■ 0, F — > cxd, with 6T ^ 1 fixed, and large temperatures, we have 



Ai^/i^p(-) ~ r(2 - 3(5T) - -T In -^ 

2 (5"(T 



(5.19) 



Since the renormalization procedure has to be performed as long as AF^'-^p^ ) < 0, we obtain a 
transition at the temperature 



T. = 



3d 



(5.20) 



where S is the parameter characterizing the distribution of quenched charges (1.6). For T < Tc. 
the renormalization procedure has to be stopped at scale T^q {T) given by 



, , STlnJ- 



3Tln J- 



2(2-3(5T) 4(1- f) 



(5.21) 



At the temperature T = Tc, the renormalization scale Teq{T) diverges : this corresponds to the 
delocalization transition found previously in References lo @ B ^3 ■ 
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C. Comparison with previous results for tiie critical temperature 

1. Case of Gaussian distribution for the quenched charges 



For the case where the initial distribution of charges is Gaussian (L7), we get 



Tc = l,=^ , (5.22) 

in agreement with the scaling obtained previously [0 |^ and with the bounds established in 
Reference [p| , which read in terms of our notations 

2(7 2cr , , 

^^- < T, < — . 5.23 

Vtt In 2 go 

We note moreover that the value of the transition temperature obtained here coincides with the 
critical temperature obtained by Stepanow et al. Iq] and is lower than the critical temperature 
obtained by Trovato et al. M], since their results read respectively in our notations : 

rriStcpanow ct al. /r 0A\ 

3qo 

/Ti Trovato et al. /r o^'^ 

'' 90 ' 



2. Case of binary distribution for the quenched charges 
For the case where the initial distribution of charges is a binary distribution, we obtain that Tc 



is the solution of the equation ( |l.ll 



ce-^«+ + (l-c)e+^«- = 1 . (5.26) 

In particular for the case q^ = q^ = qi (1.12|), we obtain 



Tc = y , , (5.27) 

which is consistent with the bounds given in Reference [pi reading in our notations : 

Tc < -4^ ■ (5.28) 

We may also compare with the bounds obtained in Reference |ll|] concerning the binary case with 
c=h,q.^ = l + h and g_ = 1 — /i where < h < 1 for which we have : 

qo^h , 2cr = l , e-2*''cosh5= 1 . (5.29) 

Our renormalization approach yields thus for this case the following behavior for the critical line 
hc{T) at high temperature 

/I \ vr /2\1 1 1 

/iJT)= — lncosh(5 L_ 2 = — ^Incosh = - + 0(— ) , (5.30) 

which satisfies the bounds established in Reference ||ll|] 

< lim (ThciT)) < 1 . (5.31) 

T — >-c>o 
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D. Thermodynamic quantities 



Following the computations done for the symmetric case in Section IVC, we obtain that the 
energy per monomer e(T) of the chain in the thermodynamic limit can be obtained as 



e(r) 



JQ.i *3i^re,(T)(*3' '') + ^re5(T)('3' '■)\ 



IqAp,^ 



/"L-TegCT) 



(Q,0 + Pr"e,(T)(Q'0] 



(5.32) 



with TeqiT) defined in Equation ( ^.2f[ ) 



Using the fixed point solution of Equations (5.1), we have at large temperature 



e{T)c^ -2(5crcoth[(5req(r)] - 
~ -2Sa - 8ad^T^y{T)e- 



sinh'[5reg(r)] 



^25r,,(T) 



+ 



(5.33) 



Since in the limit of small dissymmetry (5 — > 0, we have the relation (1.15), the energy per monomer 
at the transition Tc is simply 



e{Te) 



-qo 



(5.34) 



and coincides of course with the energy per monomer when the heteropolymer is delocalized in the 
solvent (+) for T > Tc- Surprisingly however, we find that the critical behavior near Tc is governed 
by the following essential singularity : 



e(T) - e{Tc)~ -8<j6^Tcg{T)e-^'^^^^^^ 



T-^TZ 



In^^ 

On ?!I 



(1 



T_^ exp 



In^ 

2L' 



(1 



(5.35) 



The entropy per monomer s{T) is given in the thermodynamic limit by the generalization of 
Equation (4.15) 



s{T)^ 



/q JHuM + luK - I lnZ][F+^(y)(g, + Pr,,(T)(Q' 0] 
/Q,n^rt,(T)(Q.0VPi7^(^)(Q,0f 
In l^ _ In K 



3 In/ 

= In /i = — = 

' 2 1++1 



1+ + 1- 



(5.36) 



where the averages have to be computed with the probability distributions Pp (Q, I) given in 
Equations (pTl|). Considering as before the regime S —^ 0, F ^ oo, with 7 = (5F ^ 1 fixed where 
we can use Equation ( |3.14 ), we get 



s(r)~ In/i - 12a(53Fe,(T)e-2*r,,(T) _^ _ 



3o2 In 42 
T "t- ^ 2a(l-f) 






In^ 

go 



We finally give the expression of the free-energy per monomer 

/(T)~ f{Tc) ~ 8cr(52Feg(r)e~2«re,(T) _^ i2cr53yp^^(y>)g-25re,(T)^ 



(5.37) 



~ /(re)-29o(ln— )exp 

T^Te" V qa / 



In^^ 

go 

(-ft) 



(5.38) 



which again presents an essential singularity near Tc : the delocalization transition is thus found 
to be of infinite order. 

The well known pure models presenting infinite order transitions are the two-dimensional XY 
model |2^ and the ID Ising or Potts model with inverse square long range interactions |2^ ||29|| . 
The important property of these systems is that the interaction between two defects (two vortices 
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in the XY model or two kinks in the spin chain) is logarithmic at larges distances. The usual 
estimation of the free-energy for a pair of defects, where both the energy and the entropy behave 
as the logarithm of the size L of the system, leads to a factorization of In L as in the Imry-Ma 
argument leading to Equation ( 2.15D for the heteropolymer problem. It would be interesting to 
discuss in more details the similarities/differences between these systems, but this goes beyond the 
scope of the present paper. 



E. Statistical properties of blobs 

In the localized phase T < Tc, we find that the blobs in the domain z > and z < have 



respectively the typical lengths (5.10) 



1 






1 



,^sr^.iT)_2ST,g{T)-l) 



{26T,g{T)-l + e 



-2AT.,(T)n 



As T approaches T^ , l^i^i^iT) thus diverges with an essential singularity as 



itaT) 



1 \ d-'^) 



T^T^ 4(5^(7 \6^a ^ 
whereas Ihiobi"^) ^^^^ diverges but only algebraically as 



ln-p3- 1 



a 

-2 exp 

% 


In^ 

1 ''° 


L(-i)J 


lly as 




a 


In^ 

90 



(5.39) 
(5.40) 

(5.41) 



t^t;: 



45V (1 - £) ci (i _ T ) 



(5.42) 



We note that these behaviors are compatible with the relation ( ^.13 ) between the typical lengths 
Kiob ^^'^ ^biob obtained in Reference B via an Imry-Ma argument, since we have obtained 



blob — „2 "^'6;o6 
% 



(5.43) 



Our results (5.38,5.41) also satisfy the relation ( 2.15| ) between the free-energy and the typical blob 
length Z|^gjj in the preferred solvent obtained again in Reference H via an Imry-Ma argument, since 
we have 



/(T)-/(T,)^-|(r,-r)i^ 

^ %lob 



(5.44) 



As a direct consequence of Equations ( 5.41 , 5.42 ), we obtain that the typical distances z^{T) to 
the interface in the domain z > and z < diverge respectively as T — > T~ as 



z+{T) 



Dl 



-blob 



(T) 



t^t: 



(T) - JDl-bi°\T) 



t^t;: 



/Da 

qo 



■ exp 



In 



4<T 



2(l-f 



In^ 

90 
1-^ 



1/2 



(5.45) 



(5.46) 



As in the symmetric case, the lengths l^ of the blobs are independent random variables dis- 



tributed respectively with the probability distributions Pp ,j^\{l)- In particular, near Tc, the 



rescaled variable 



A, 



4,^2^- 



l 



,2Sr,g{T) 



T^T- 





1+ 


t^-p 


'"ff 


W-m 



(5.47) 



is distributed with a simple exponential law e + (see Equation 5.14) 
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F. Density profiles p (z) 

Following the calculations done previously in the symmetric case, we obtain that the probabilities 
p^{z) that a given point of the polymer is at a distance z from the interface in the (±) domain 
(with the normalization J^ dz{p'^{z) + p^{z)) = 1) reads 



p±(z) 



!^dii{p^{i) + p-{i))j. 



dlP^{l) 4 




due 



(5.48) 



In the regime ^F » 1 studied above where P^ Q-) has the asymptotic simple form of Equation 
(5.14), we get the following scaling form 



p+{z)dz = R+{Z)dZ 
where the rescaled distance Z to the interface and the scaling function R{Z) read 

iSz z 

Z = 



/2X)g5r,,(T) j,^j,- r^ 



2 go 



exp 



In 42. 



^^W) 



(5.49) 



(5.50) 



and 



i?+(Z) =4 



due 



dvv e 



2„-V 



(5.51) 



Deriving with respect to Z yields 
dR+{Z) 



dZ 



-4:Z' 



due^'^^^^X; = -AZKi{2Z) 
u-^ 



(5.52) 



In particular, the asymptotic behaviors at Z ^ and Z ^ oo are given respectively by 



R+{Z)^--2Z + 0{Z^) 



(5.53) 



and 



R+{Z) ~ VnZe 

Z—*-oo 



-2Z 



(5.54) 



G. Correlation function < sgn(2;(s))sgn(2;(s')) > 



Within the renormalization picture, the correlation function < sgn(2;(s))sgn(2;(s')) > of the sol- 
vents seen by two monomers at distance {s — s') along the polymer chain containing a dissymmetric 
distribution in hydrophilic/hydrophobic components corresponds exactly to the spin correlation 
function {So{t)Sx{t)) of the random field Ising model in the Glauber dynamics starting from a 
random initial condition in the presence of an external field, which is computed in Reference pit] . 
Here we do not give the full result, but only the correlation length ^(T) which, not surprisingly, is 
simply given in the regime 6T ^ 1 that we consider by the typical behavior of l^iobC^) obtained in 
Equation KM) 



aT)^lL,iT) 



-TT exp 



In^^ 

90 



1-^ 



(5.55) 
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VI. FINITE SIZE PROPERTIES 

Up to now we have always considered the thermodynamic hmit of heteropolymers of infinite 
length L — > oo. In this section, we study the finite-size properties of the localization at the interface 
for the case of periodic boundary conditions, i.e. we consider cyclic chains (gi, (72, ■ • ■ 9l, 9l+i — Qi) 
of finite (but large) size L, with positions (zi, . . . , zl+i = zi). We can apply the renormalization 
procedure as before, but the procedure will now stop at some finite renormalization scale rdcioc. 
where the chain delocalizes from the interface. This renormalization scale depends upon the 
realization of the quenched charges and corresponds to some temperature T^eioc via the relation 



dcloc 



reg(Tdcloc) (6.1) 



given in Equations (4.10) and (5.21) respectively for the symmetric case and the biased case. In 
the following, we characterize the distribution of Tddoc for the ensemble of finite chains of size L, 
and the probability distribution of the number of blobs for the chains of length L that are still 
localized at a given temperature. 



A. Probability measure for the blobs of a cyclic chain of length L 

As shown in References |lq|, [M, it is possible to follow the renormalization procedure for 
finite size systems. Here we deal with periodic boundary conditions, and to avoid problems with 
the translation invariance along the chain, it is convenient to mark a point of the chain, called 
"the origin" in the following, and to give an orientation to the chain. For fc = 1, 2 . . ., we define 
-^r L (Qij ^1! Q2, Wt ■ ■ Q2k, hk) as the probability that in the chain of length L at renormalization 
scale r, the origin belongs to a blob (Qi, li) in the domain (±), and that there are exactly (2fc— 1) 
other blobs in the chain of absolute charges and lengths given by the sequence {Q2, h', ■ ■ ■ Qik, hk)- 
We also need to introduce the probability Nj.\ (Q) that the chain at scale F is already delocalized 
in the (±) solvent with absolute charge Q. The normalization of these probabilities read 

V / N^''+{Qi,h;Q2,l2;...Q2k,l2k) 

00 „ 

+ E / N^''LiQlJl-,Q2,l2;...Q2k,l2k) 

+ / dQNl\{Q) + / dQN^-^iQ) = 1 . (6.2) 



Note that in the last two terms, the absolute charge Q is a random variable of the domain [0, +00 [, 
contrary to the other terms where by definition of the renormalization rule we have Q,; > F. 
The renormalization equations for these probabilities read for fc = 1, 2, . . . 

^^^r,L {Qi^h',Q2-,h', ■ ■ ■Q2k-,hk) = 
2k 

y. / -^r L ' (Qi: ^1; • ■ ■ ; Qj-ii ^i-i; Qj ^; r, ' ; Q J ',Qi+i,li+i---',Q2k:l2k) 

i=2 J Q+Q" -'^=Qi,i+i' +1" =h 

^T L ' (Qj^S-T, ^;(5 ,1 ',Q2,h---',Q2k,hk) 
Q+Q"-r=Qi,l+l'+l"=li 

jQ+Q"-r=Qi,l+l'+l"=h 

f N^''+''^{r,l';Q"J";Q2j2;...;Q2kJ2k:Q^l) > (6-3) 

Jo+o"-r=o^.i+i'+i"=h 



and 






h) ■ (6.4) 
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As already obtained in References |g5[ [|9[ for the case of fixed boundary conditions, the above 
renormalization equations concerning periodic boundary conditions are solved by a quasi-factorized 
form for k — 1,2,... 



^T.L {Qlih]Q2,h',- ■ ■Q2k,hk) 



2k 



= iiP^{QiM)P^{Q2M) ■ ■ ■ PT{Q2k-i.i2k^i)P?{Q2kMk)5 \L-Y,h 



(6.5) 



where Pp (Q,/) are the bulk distributions given in Equation (O), satisfying 



drP^{Q,l)= f 
JQi 



+Q3-T=Q.li+l2+la=l 



P^{T, l2)P^{Ql,ll)P^{Q3, h) 



/•OO 

+p^{Q,i) / di'{p^{T,r)-p^{r,n) 

Jo 



(6.6) 



Note that in the measure (p]9), the first bond plays a special role as it is defined as the bond 
containing the origin. 



The equation (6.4) now reads 

drN^iiQ)^ f hP^{Q + r,h)P^{r,l2)S{L-{h+l2)) 

Jll,l2 

+ f hP^iT,h)P^iQ + r,l2)5iL-ih + l2)) 



= LPr±(Q + r,.)*LPrT(r,.) 



Using now the explicit expressions (5.1), we get in Laplace transform with respect to L 



'dLe-P^(9riV^;±(Q)) = -dp (- 



sinh^[rA(p)] 



~Q(A{p)cothlrA{p)]TS) \ 



where A(p) = ^S'^ + f . 



B. Probability for a chain to be delocalized at renormalization scale F 



(6.7) 



(6.8) 



The equations ( |6.7D for the probabilities N-p'i^ (Q) that a chain of length L is already delocalized 
at scale F with an absolute charge Q in the domain (±) can be integrated with respect to F to 
yield 



Jo '^r'^^^^-2alv A(;,) sinh^[rA(p)] j ' 



(6.9) 



rl,±/ 



We may check that A'p'^(Q) vanishes as it should for F ^ 0. In the limit F -^ oo, the chain is 

expected to be always delocalized, with a charge corresponding to the sum X]i=i 1i ^^ ^ independent 
identical random variables. Indeed, using ( 1.15| ), we find 



Nrt^.L{Q) = c;i^ 



-Qiy/S^TS) 



2syF+f 



1 (QT-C-gp)^ 



2V2ctL 



(6.10) 



rl,± 



as expected from the Central Limit Theorem. 

For arbitrary F, we can obtain the total probabilities N-^'^ that the chain is already delocalized 
at scale F in the domain (±) as 






'-p^L 



dQN^ltiQ) 

1 /cosh[FA(p)] F \ 

2ct(A(p) cosh [FA(p)] t S sinh [FA(p)]) 1^ A{pj sinh[FA(p)]y 



(6.11) 
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In particular for the symmetric case 5 = 0, we have 



N, 



i,± 



aL\ 



r..-2AA^A.-l 



(6.12) 



The scahng function Af{X) represents the total probability for a symmetric chain of length L to 
be delocali zed a t temperature T corresponding to the renormalization scale T = T^q (T) defined in 
Equation (4.1C). It reads 



AA(A)= c:i, Q 



3inh(2Vs) 



+ 00 



EM) 



'e-^" "t- =1- 



h^ 



+ 00 

E 



(6.13) 



1. Distribution of the delocalization temperature for symmetric finite chains 

The probability that the polymer chain delocalizes from the interface between T and F + dF is 
given by {dr [Af{X)] ) . To characterize the distribution of the delocalization temperature Tdcioc for 
the ensemble of the symmetric finite chains of length L, it is thus convenient to define the rescaled 
variable 



^ eg\-^ dclocj 



(6.14) 



where the function Feg(T) has been defined in Equation (4.10). The final result is that g is 
distributed with the law : 



2 +00 



^n e ^s'-' 



= ^ Y. [2(2m + l)V-l]e-(2'"+i)'s' 



(6.15) 
(6.16) 



2. Distribution of the delocalization temperature for dissymmetric finite chains 

In the biased case, the probability that a chain of length L delocalizes from the interface between 
F and F + dF in the domain (±) is given by 



drN^i = Cr^ 



■p^L 



-d„ 



A2(p) 



^ Vsinh[FA(j5)](A(p) cosh [FA(p)] T (5sinh[FA(j5)]) 



(6.17) 



In the regime 7 = (5F ^ 1 considered before (see Equation (5.14)), the inverse Laplace transform 
is dominated by a single pole contribution 



drNp+ ~ 2(5a+(7)Le-^"^(^) 



(6.18) 



where a~^{'~f) = 4(T(5^exp(— 27) has been introduced in (5.14). We now use the correspondence 
between temperature and renormalization scale given by F = Feg(T) (5.21). To characterize the 



distribution of the delocalization temperature Tdeioc for the ensemble of dissymmetric finite chains 
of large length L, it is thus convenient to define the random variable 



1 



1 



46^ aL \d^<T^ 
The final result is that r is distributed with the law 



(6.19) 
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D+(r) = le-^ . (6.20) 

In particular, the typical value for the delocalization temperature corresponds to r ~ 1 and thus 
presents a correction in l/(lni) with respect to the critical temperature Tc 



«c-Tc(l-^) . (6.21) 



C. Distribution of the number of blobs at scale F 

The probability that a chain is still localized around the interface at scale F and is in a state 
with (2fc) blobs read for fc = 1, 2, . . . 

Bf{T) = / Nl\+{Q^M;Q2, h;.-- Q2k, hk) 

+ [ Nl\~{QiM;Q2M-:---Q2kMk) ■ (6.22) 

The Laplace transform of the generating function of the probabilities B^ (F) reads 



f "i^-"' g.^Br(r,) ^ [-oAP-MP^iM j-^^ 



(p) 



Zi 



sinh[rA(p)] ((J^cr sinh[FA(p)] +pFA(p) cosh[FA(j5)]) 



(6.23) 



(52(7(1 - z) +p(cosh^[rA(p)] - z))(,5V+pcosh^[rA(p)]) 
In particular, the average number of blobs reads 

^^-1 /^ <?V + prA(p)coth[rA(p)] \ 

. J£i^ + £ ,-(.-.*^.4_£^?!4_ f 1 + ^^^L- - 2^) . (..24) 

sinh^((5F) ^^ F2(n27r2 + (52F2) ^ n27r2 + (52F2 F2 / ^ ' 



The first term is proportional to L and corresponds of course to the ratio —^— that could have been 

anticipated from the study of Section M concerning the thermodynamic limit i ^ oo. The other 
terms, which represent the finite-size corrections to this dominant contribution, decay exponentially 
with L. 

We may also compute the second moment 

< (2fc)2 >,,p^ £4^252. (F) = C~l^ U [^d, (Pp+(p)Pp-(p))] -' + ^r+(P)^r- (P) 



= c-L,U 



il-P^ip)Pfip)f 

(2,52cr+p(cosh^[FA(p)] + l))((52g+pFA(p)coth[FA(p)]) \ 

"P^^ V" p3sinh*[FA(p)] J 

The dominant behavior of the variance of the number of blobs is thus given by 

< (2fc)2 >L,r -(< 2fc >L,rf = -^ (cosh(2(5F) - A5T coih{5T) + 3) + ... (6.26) 

sinh (5F) 
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In the symmetric case 5 = 0, the generating function may be more expHcitly computed as 
f> k„2k(j.,r-i ( rztanh[rA(p)] \ 

h '^' V ^^^^'^" [rAb)] - z) ) 

+00 

= E (e^7*("+"-)'-e-7*(f +"-)') , (6.27) 



n— — 00 



where a = ArcCos-yz e (0, f ) for z G (0, 1). In particular, the average and the variance of the 
number of bfobs read : 

< 2fc >,.,r= 2^ + E V ^ " V" '^ e 7^" '^ (6.28) 

n=l ^ ^ 

< (2A;)2 >i,r -(< 2fc >i,r)2 = ^ + ... , (6.29) 



2^ 



where F = reg(T), and where the dots represent terms that are exponentially small in {L/T 



VII. NON-EQUILIBRIUM DYNAMICS STARTING FROM A ZERO-TEMPERATURE 

INITIAL CONDITION 



As explained in Section IIIB, the dynamics at temperature T for i > starting from a zero- 
temperature initial condition at t = 0, i.e. after a quench to T = for i < 0, can be described by 
the renormalization procedure, where the renormalization scale F now corresponds to time via 

F(i)=rin-^ (7.1) 

to 

as in References |lq| p9[ . The dynamics takes place up to time t^q where equilibrium at temperature 
T is reached : 

Tln^=re,(r) . (7.2) 

to 

In the symmetric case, we thus have the scaling 

teq^Tl^{T)^T\\nTf , (7.3) 

whereas in the biased case 

/ 31n4- , 
ie, ~ exp ^^, . (7.4) 

All the quantities computed before for the equilibrium at temperature T as functions of Fg^ (T) 
have the same expressions for the dynamics as functions of F(i) = Tlni for large time t < t^q- For 
instance, in the symmetric case, the typical length of blobs behaves as 

a a 

whereas in the biased case, the typical lengths of blobs in the domains z > and z < behave 
respectively as 






4,5V 4(52ct -2' 



£W Tint _ Tint 
^'^ 25a ^ 25a " -70 ' ^ ' 
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VIII. CONCLUSION 

In this paper, we have proposed a new approach based on a disorder-dependent renormahza- 
tion procedure to study the locahzation of random heteropolymers at the interface between two 
selective solvents within the model of Garel et al. |g]. The renormalization procedure has been 
defined to construct an effective thermodynamics, where one only considers the heteropolymer 
configurations that correspond to the optimal Imry-Ma domain structure. At high temperatures, 
where the distribution of absolute charges of blobs generated by the renormalization procedure 
becomes infinitely broad, the effective thermodynamics is expected to become accurate and to 
give asymptotic exact results. With the renormalization approach, we have recovered that a 
chain with a symmetric distribution in hydrophilic/hydrophobic components is localized at the 
interface at any finite temperature in the thermodynamic limit, whereas a dissymmetry in hy- 
drophilic/hydrophobic components leads to a delocalization phase transition, in agreement with 
previous studies. In addition, for both cases, we have given explicit expressions for the high tem- 
perature behaviors of various physical quantities characterizing the localized phase, including in 
particular the free-energy per monomer, the distribution of the blob lengths in each solvent and 
the polymer density in the direction perpendicular to the interface. For the case of a small dis- 
symmetry in hydrophobic/hydrophilic components, where the delocalization transition takes place 
at high temperature, the renormalization approach allows to study the critical behaviors near the 
transition : we have found that the free energy presents an essential singularity at the transition, 
making all of its derivatives continuous at Tc (infinite order transition), that the typical length of 
blobs in the preferred solvent diverges with an essential singularity, whereas the typical length of 
blobs in the other solvent diverges algebraically. We have then studied the finite-size properties of 
the problem by considering cyclic finite (large) chains. In particular, we have given the probabil- 
ity distribution of the delocalization temperature for the ensemble of chains of fixed finite length 
L, and the distribution of the numbers of blobs in the chain still localized at some temperature. 
Finally, we have briefly discussed the non-equilibrium dynamics at temperature T starting from a 
zero-temperature initial condition. 

In conclusion, the disorder-dependent renormalization approach yields explicit predictions for the 
high temperature behaviors of various physical quantities, that would certainly be interesting to test 
precisely by numerical studies. In particular, the predictions concerning the critical behaviors near 
the delocalization transition are certainly surprising, and disagree with the Monte Carlo studies of 
Reference B where power laws were found, and with the Replica Gaussian variational approach of 
Reference M] where the transition was found to be of second order. It would thus be particularly 
worthy to test numerically more precisely the order of the transition. We hope that the predictions 
of the renormalization approach concerning the finite-size properties of the problem presented in 



Section VI will make easier the comparison with numerical studies. Finally, we would like to stress 
that the renormalization approach yields not only predictions for averages over the quenched 
disorder, but also precise predictions sample by sample. For instance, for a given realization 
of the quenched charges, one may compare the typical configurations of the heteropolymer at 
equilibrium at temperature T with the blobs structure obtained by the numerical implementation 
of the renormalization procedure up to scale Teq{T). For the case of cyclic chains, one may also test 
numerically sample by sample the renormalization prediction for the delocalization temperature. 
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APPENDIX A: IMRY-MA ARGUMENT FOR THE PROBABILITY DISTRIBUTION 

OF BLOB LENGTHS 

In this Appendix, we recall the determination via a direct Imry-Ma argument ||25|1 of the asymp- 
totic behavior of the probability distribution of blob lengths in the limit of small lengths, and 
compare it with the renormalization result. 
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As noted by Garel and Orland pSl , one may extend the usual Imry-Ma argument for the het- 
eropolymer exposed in Reference [^ to obtain not only the scaling of the typical length of Imry-Ma 
domains, but also some information on the probability distribution of the length of domains. The 
idea is as follows [^. On one hand, the energy associated to a blob of length I is the sum of / 
independent random variables and is thus typically given by 

El = -uV2^l , (Al) 

where u is a random variable of order 1 , that has at least to be positive since the system tries to 
use the favorable fluctuations of the disorder to lower its energy. According to p5[ , the probability 
distribution for large u is simply given by the Gaussian of the Central Limit Theorem : 

G{u) - e-"'/2 . (A2) 



^oo 



The behavior of the distribution G{u) for small u seems however not simple to determine a priori, 
because the small values of u are certainly suppressed with respect to the Gaussian distribution. 
Indeed, if an Imry-Ma domain corresponds to a small value of u, it is very likely that the het- 
eropolymer will prefer an organization into other Imry-Ma domains. It seems however difficult 
to write precisely the consequences of this effect on the distribution G{u). On the other hand, 
the entropy associated to a blob of length I behaves as (— |ln/) (in Q, only powers of / were 
considered for the symmetric case and thus In I was replaced by ^" ~ 1, but here for a better 
comparison with the renormalization approach, it is more convenient to keep this In / dependence) . 
The minimization of the free energy per monomer 

m^-J^^+l^'^ (A3) 



with respect to I leads to 



ST Inl 



(A4) 



(A5) 



V2^ Vl 
The typical blob length ^*^p corresponding to utyp ~ 1 thus satisfies 

VI*w 3T 

InitvP^ ^/2^utyp 

p2 //TIN 

Since within the RG approach 1*^^ = -^ — , the above equation corresponds to equation (4.10) 
with the choice Utyp ~ v2 i.e. a numerical factor of order 1. Now introducing the rescaled length 
of Imry-Ma domains 

^ = ii ' (A6) 

we find that it is simply related to the random variable u through 

.^-L . (A7, 

The behavior (A2) of the distribution G{u) at large u thus corresponds to the following behavior 



of the distribution P{X) at small A : 

PW-o^e~^ (A8) 

where C is a constant of order 1. This asymptotic behavior corresponds exactly to the behavior 



(4.2C) found via the renormalization approach. Concerning the behavior of G{u) for small u, we 



explained above why it is difficult to determine it a priori. Using the relation (A7), the asymptotic 
behavior ( [4. 191) found via the RG approach 

PrgW- e-^^ (A9) 

A — *oo 

would correspond to 

G{u) - 4e"^ (AlO) 

i.e. a very strong suppression of small values of u for the energies of Imry-Ma domains. 
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